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This paper is concerned with the non-linear vibration performance of a rigid 
rotor supported on three bearings, one being surrounded by a squeeze -film damper. 
This damper relies on the pressure built up in the squeeze film to help counter- 
act external forces arising from unbalance and other effects. As a result a 
vibration orbit of a certain magnitude results. 

Such vibration orbits illustrate features found in other non-linear systems, in 
particular sub-harmonic resonances and jump phenomena. Comparisons between theoreti- 
cal predictions and experimental observations of these phenomena are made. 

INTRODUCTION 

The rotors of aero-engine gas turbines are often supported on three rolling- 
element bearings, the intermediate bearing usually providing a thrust capacity. The 
reason for this is that axial expansions take place from the thrust face and, being 
proportional to distances along the shaft to the ends of the rotor, these expansions 
do not impose such stringent demands on axial clearance at these ends as they would 
if the thrust bearing were at one extreme end of the rotor. 

However, many aero-engine rotors are stiff, (running well below any flexural 
critical speed), and this means that problems of misalignment almost inevitably 
exist. The presence of a squeeze -film damper around one of the bearings helps 
considerably to alleviate these problems. However, since such dampers carry very 
little radial load under these circumstances, non-linear response to unbalance and 
sub-harmonic resonance with the casing flexibility are likely to occurs 


Non-linear response manifests itself as persistence of synchronous vibration 
well beyond the undamped natural frequency of the assembly with a sudden jump-down 
to a lower vibration level at some higher speed. During engine run-down the lower 
vibration level persists until a speed lower than the original jump-down speed is 
reached, when a sudden jump-up will occur in the vibration amplitude. On occasion a 
jump-up with run-up has also been observed. 

Sub -harmonic resonances occur due to the weak non-linearity offered by the 
squeeze-film damper, usually at rotational speeds equal to some integer multiple of 
the natural frequency of either of the bounce modes of the rigid rotor in its 
casing. These resonances can be just as damaging in their own way as the resonances 
which occur at speeds equal to these natural frequencies (that is when the integer 
is unity). 
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FEATURES OF THE PRESENT INVESTIGATION 


To elucidate these problems, consider a stiff rotor in its flexible casing 
(Fig. la), in which the intermediate bearing is surrounded by a squeeze- film damper. 
To formulate ideas let the stiffnesses of bearing pedestals 1 and 2 be high compared 
with that of bearing 3. Further, assume that bearing 2 is held fairly centrally 
within its squeeze -film damper. Under these circumstances the equivalent system can 
be drawn as in Fig. lb, in which m, k and b are equivalent properties decided 
by the dimensions of the original assembly. in particular, b will describe a 
non-linear damper. 

During running any vibration is likely to produce a concentric orbit and the 
squeeze-film damper is not called upon to provide a static load -carrying capacity. 
Consequently, if the oil -supply pressure is of reasonable magnitude, the damper 
annulus remains full. The damper is often supplied with oil to a central circumfer- 
ential oil groove which effectively separates the squeeze film axially into two 
parts or lands. Governing parameters are usually assumed to be appropriate to one 
land. 


If the oil supply to the damper is sufficient for the squeeze -film to remain 
full and uncavitated a rotor vibration response curve such as (a) in Fig. 2 may be 
expected [1], [2]. If however cavitation occurs, this has the effect of introducing 
a stiffness term in the radial restoring force provided by the squeeze film and 
curves such as (b) are predicted [1]. The latter curve reveals a hardening - 
stiffness non-linearity and a deleterious influence on rotor response. At moderate 
frequencies three amplitude ratios are possible, the intermediate ratio being 
unstable, resulting in a jump-down with engine run-up and a jump-up with engine 
run-down. Such jumps have been reported in practical systems, as has the occasional 
jump-up with increase in engine speed. Fig. 2c shows the result of accelerating a 
rotor up to a speed of 3500 rev/min. and back over a period of 8.3 minutes. The 
degree of blackness depicts the strength of each harmonic from 1 x rotor speed to 3 
x rotor speed, the other harmonics being practically non-existent. A predominant 
feature is the disappearance of the 2 x and 3 x harmonics after the jump-down on 
run-up and a reappearance of them after the jump-up on run-down. This suggests that 
the jumps form a demarcation between effectively linear and non-linear vibration and 
have been noted in test engines. Numerical investigations have also been carried 
out to test the possibility of sub-harmonic resonance arising [1]. Fig. 3 shows 
predicted sub-harmonic transients for a rotor-bearing assembly incorporating a 
squeeze-film offering fairly light damping. Subharmonic vibration occurs at speeds 
in the regions of two, three and four times the lower bounce mode critical speed. 
Such effects have also been reported in practical systems [3]. 

Very few experimental programs have been conducted to investigate problems of 
non-linear jumps and subharmonic resonances in rotor -bearing assemblies. The experi- 
ments of Nikolajsen et al [3] and Simandiri et al [4] are limited respectively to 
flexible transmission shafts with squeeze-film isolators and to a highly idealised 
test facility, showing little resemblance to an engine assembly. 

RESPONSE CALCULATIONS 

Fig. 4 shows in diagrammatic form the outer race of a rolling -element bearing 
within the oil container (the bearing pedestal), under the action of a restoring 
force kce arising from the equivalent stiffness k of a parallel spring as in 
Fig. lb. Vibration results from a dynamic force P c due to some cause, such as 
unbalance. The amplitude of orbital motion depends upon kce, P c , p-j and P 2 , the 
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last two forces arising hydrodynamical ly from the squeeze-film. Any gravity load on 
the squeeze-film is assumed to be effectively neutralised by the parallel spring. 


The equations governing the concentric motion of the shaft centre are then, 
P c cos (wt — vp ) — Pi - k.ce = — mcew 2 

« • I 

P c sin ( cot — \p ) — P2 = 0 


• • • ( 1 ) 


where e is the dynamic eccentricity ratio resulting from unbalance. In equations 
(1) the squeeze-film forces, Pi and P2 , can be shown from hydrodynamic considera- 
tions, to be given, for concentric motion, by 


Pi = 0 

P2 = - 


TTT)Rfi3 

c 2 ( 1 — e 2 ) 3/2 


Equations ( 1 ) may be made non-dimensional by dividing by mew 2 to give 
Q c cos (wt— vj/) = — e(l — R) 

and 

~ . n-r|R ,8.0 e 

Q c sin (wt - *) - — (-) 3 - (1 ._ e Ji) 3/ 2 

where K = k/mw 2 = (w/w n )~ 2 , and Q c = P c /mcw 2 . 

If P c arises due to mass unbalance mu, then P c = muw 2 and Q c = u/c. 
TiRf fl/c 

Putting (3 = — — , equations ( 3 ) may be re-written to give 

mo) 

— cos (wt— vp ) = ( K — l)e 


— sin (wt — 4/) = 7 T/ 3 e/(l— e 2 ) 3/2 


Hence, after some manipulation, we obtain 


where 


-l 2 

(1- e 2 )3 

f- r = a2 

f - ] 2 + [i - 

f “I 

c. 


1 fc> n j 

l “n J 

l “n J 


w n = ^k/m 

Assuming a fairly typical value of A = O.O 7 , say, a family of response curves 
of versus w/w n may be constructed for different values of u/c (Fig. 5 a). 

From these curves an indication can be obtained of the amount of vibration e, 
suffered by the rotor in the symmetric bounce mode, a common mode of vibration in 
many rotating assemblies. Of equal importance is the dynamic force F transmitted to 
the engine frame. Bearing in mind that Pi is zero, the force F is given by 

F = [P2 2 + (kce) 2 ] 1/2 ... (6) 
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and is best expressed as a ratio (the transmissibility) of the unbalance force P c . 
This transmissibility T is given by 


T = F/P c = F/mctj2 
u/c 

From equations (2), ( 6 ) and (7) we obtain 


T2 = [ (i-^T 2 )3 + k 2e2 ] /(«/c)2 


Now from equation (4), 

nZp2 € 2 

(1 - e2 )3 


= (u/c )2 


e 2 [k — i] 2 


(7) 


and so finally 


T 2 = 1 + (— — )2 

v u/c ' 


[ ( w /« n )2 11 


Again, assuming a value of A of 0 . 07 , a family of curves showing the dependence of T 
on u>/w n for different values of u/c is shown in Fig. 5b. 


An inspection of Figs. 5a and 5b shows that both e/(u/c) and the transmissibi- 
lity T depend upon u/c. This is ejqpected since the system is essentially non-linear. 
Also, only for w/u> n above ^ 2 , is T less than unity. This speed ratio should thus be 
made to correspond to the lower limit of normal operation of the engine, when it can 
readily be seen from Figs 3a and 3b that a low value of u/c is desirable on three 
counts - low vibration amplitude, low transmissibility and low transmitted force. 


Now suppose that the supply pressure is insufficient to maintain a positive 
pressure in areas of the squeeze film where the boundary surfaces are instantane- 
ously separating. For the sake of simplicity assume that as a result the squeeze 
film becomes half cavitated. It may then be shown that P 2 is halved and that Pi is 
no longer zero, being given by 


_ 2r)Rl3w e 2 

Pl ' • (i---i7)2 

Hence, for circular concentric whirl, we have the following equations which 
correspond to equations (4) for the uncavitated case 


u — 2 j 3 e^ 

— cos (u>t— v|>) = (k - l)e + (1 e2)2 

u 7773 e 

- sin(wt— Y) - . (1 _ e 2)3/2 


(8) 


- \ 

in which /3 = A( k ) /tt . 
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+ ZK - 1] 


Using equation (9) this may be rewritten as 


in which R 


= k/mo)2 = 

0) n 


-2 

) 


The effect of cavitation is shown in Figs. 5c, d, which reveal a deleterious 
influence on vibration and on transmissibility . The distortions in the response 
curves of Fig. 5c compared with what one would expect from a linear system are 
responsible for the so-called jump phenomenon. At high frequencies three amplitude 
ratios are possible, the intermediate one being unstable. 

STABILITY CONSIDERATIONS 

In most applications of a journal bearing we are concerned with the stability 

of a static equilibrium position e.g. [6]. In the present work it is necessary to 

consider the stability of a circular concentric equilibrium orbit . The two cases can 

be compared as shown in Fig. 6a, b. It may be shown [6] that "for perturbations from 

a static equilibrium position the hydrodynamic forces P^ and P 2 of a journal bearing 

oil film are given by 

/ . S/Wi 

Pi ,2 = XWj(l— 2a/Wj )fi ,2(e, y _ > 

With little loss of accuracy [7] we may write them as 
Pi = eFi(e) ± — a) hi(e) 

p 2 = eF2<e) ± ( yj-— a) h2<e) 

in which the positive sign relates to the squeeze damper and its force directions 
(Fig. 6b) and the negative sign to the journal bearing and its force directions 
(Fig. 6a). Now in the case of the squeeze damper, if the journal centre is moving 
on a circular concentric orbit at angular velocity a s , then Wj = 0 and 5 = u> + a s . 

Hence we finally have for the journal bearing 

Pi = eF'i ( e ) — ( - a ) hi ( e ) 

P 2 = eF 2 (e)- ( — a) h 2 (E) 

while for the squeeze -film damper 

Pi = eFi ( e ) — ( a s + w ) hi ( e ) 

P2 = eF 2 ( e ) — (a s + w)h 2 (e) 

Now for a journal bearing, for a « t»)j and neglecting second order terms 

dP r = dPi — P 2 da = F^de — — de + h^da + — h 2 da 

*“ 2 OS 2 

and , 

t t * Cl) -i dh^ • CO -i 

dP s = dP 2 + Pi da = F 2 de — — 1 — 2 de + h 2 da — — hi da 

2 OS 2 
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while for a squeeze damper, for a s « « 

dP r = d Pi + P 2 da s = Fide — u ^j^de - hida 3 — <oh 2 dot s 

/ / ah • • • ( 10 ) 

dP 3 = — dP2 + Pi da s = — F2de + de + h2 da 3 - whi da 3 

Now r, s correspond to the e and a directions, respectively. Thus for the squeeze 
damper the direct damping terms are equal and of the same sign as for the journal 
bearing; the cross damping terms are equal but of opposite sign; the direct 
stiffness terms are twice the values and of the same sign, and the cross stiffness 
terms are twice the values and of opposite sign. 

The effect of adding an extra radial force kce due to the retainer spring is to 
increase dP r and dP s respectively in equations ( 10 ) by kcde and kc eda s . This in 
turn has the effect of adding the stiffness X to each of the direct stiffness 
terms in equations (10). Applying these new equations to perturbations of a rotor 
from a concentric circular orbit indicates that part of the response curve is indeed 
unstable . 


SUBHARMONIC RESONANCE 

Figure 3 shows subharmonic transients, which could persist if some negative damping 
were present to counteract the positive damping from the squeeze film damper. Such 
negative damping could emanate from any journal bearings present [3] or from inter- 
nal shaft friction, say. 


TEST RESULTS 

To create a realistic engine configuration, a two-bearing test facility used in 
a previous research programme [5] has been adapted to form a three -bearing rigid 
rotor incorporating the essential features of a medium-sized aero-engine assembly.lt 
is not uncommon in gas turbine rotor vibration for an antisymmetric rigid-rotor mode 
to exhibit a node near one of the rolling -element bearings. With such an 
application in mind, the test rig shown in Fig. 7 was used for the investigation of 
the squeeze-film damper, 1, at one of its three rol ling -elements , 2. The self- 
aligning bearing, 3, constituted a pivot about which an antisymmetric (i.e. conical) 
mode of vibration would occur, when the rotor, 4, was acted upon by a force arising 
from rotation of the unbalance mass, 5. Flexible bars, 6, simulated the pedestal 
flexibility of an actual engine. The flexible bars were screwed into a heavy 
foundation block, 8, which represented ground. The squeeze-film dimensions sere set 
by the outer diameter of the rolling bearing (136 mm), and damper land (9 mm) and 
its radial clearance (0.216 mm). Oil of 6 cp viscosity was supplied to the damper 
via three supply holes, 9, and a central circumferential groove. End plates were 
attached at the ends of the outer element of the damper to afford some sealing. 
Proximity vibration pick-ups were used to measure the vibration of the shaft 
relative to the pedestal and relative to ground, and strain gauges were provided on 
the flexible bars from which transmitted forces were recorded. 

Tests indicated the existence of jump phenomena during run-up and run-down, 
when the squeeze -film damper ran at low supply pressure (Fig. 8a), but no such jumps 
when it can with a higher supply pressure (Fig. 8b). These results agree with the 
trends predicted by Figs. 2b and 2a and by Fig 5c, d. 
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Figs. 9a, b show further comparisons, this time between the predicted and the 
measured horizontal components of force transmitted to ground, as measured by strain 
gauges fixed to the flexible bars 6 of the rig (Fig. 7). The same trends can be 
seen in both, the one exception being the continuous fairly steep rise in trans- 
mitted force with speed for the largest unbalance (1.9 gr-m) in Fig 9b* Fig. 10a 
shows a predicted waterfall diagram which illustrates vibration occurring at 
engine -orders half, one and two. These engine orders are encouraged by such factors 
as transients or slight eccentricity from the centre of the squeeze -film damper. The 
rig's natural frequency component at 33Hz is however present at all rotor speeds, 
and, in particular, is excited by the half -order when the rotor speed is 66 rev/min. 
The jump-down in the 1EO is clearly discernible and coincides with the disappearance 
of the 2EO. As such it agrees with the experimental observations illustrated in 
Fig. 2c. 


Finally, Fig. lOb shows an experimental waterfall diagram, for comparison with 
Fig. 10a. Common features can be clearly seen, such as the strong first engine 
order, the %E0 exciting the natural frequency of 33 Hz at a rotor speed of 66 Hz and 
the 2EO, with its disappearance after the jump in the 1EO. Owing to the scale 
adopted to bring out the %EO and 2EO, the 1EO was everywhere extremely large, but 
the jump certainly occurred at about 48 Hz as can be verified from Fig* 9b. The 
response at 84 Hz illustrates another natural frequency of the experimental rig 
which was not catered for in the theoretical model . 

CONCLUSIONS 

This paper has considered the case of a flexibly supported three -bear ing rigid rotor 
utilising a squeeze-film damper at one of its bearings. It has shown that the 
assembly can be regarded for analysis as a parallel combination of spring and 
non-linear damper supporting the rotor mass. As a result, subharmonic resonances 
and jump phenomena are predicted. 

A test facility has been described, which displays such phenomena and these have 
been presented in the form of response curves and a waterfall diagram. The fre- 
quency content in these shows close resemblance to theoretical predictions presented 
in the same form. 

Finally, a relationship has been established between the linear stiffness and 
damping coefficients of a rotating- journal bearing for perturbation from an equilib- 
rium position and the corresponding coefficients of a non-rotating squeeze-film 
damper for perturbation from an equilibrium orbit. 
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Figure 1. - Rotor and equivalent system. 
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(a) Uncavitated rotor vibration response curve. 

(b) Cavitated rotor vibration response curve. 



(c) Frequency diagram. 

Figure 2. - Response curves and frequency diagram. 
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Figure 3. - Subharmonic transients; A = 0.014, u/c = 0.091; w film 













(a) Journal bearing; r and s are 
stationary coordinates . 



(b) Squeeze-film damper; r and s are 
rotating coordinates. 


Figure 6. - Comparison of journal bearing 
with squeeze-film damper. 
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6 Flexible bars 

7 Housings 


1 Squeeze-film damper 

2 Rolling element bearings 

3 Self-aligning bearing 

4 Rotor % Foundation block 

3 Unbalance ma^s 3 Oil supply 

Figure 7. - Experimental rig. 



(a) Supply pressure, 2 lbf/in 2 . (b) Supply pressure, 10 lbf/in 2 . 

Figure 8. - Responses. 



1>TTTJ ! 48 rotor speed, Hz 

(b) Measured responses. 

Figure 9. - Horizontal components of force transmitted to ground. 


363 








